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A, 


^-* 1 + 

2 


Q 


(14.54) 



This is energy equation, which is valid even in the presence of friction or non¬ 
equilibrium conditions between Secs 1 and2. Itis evident that the sumof enthalpy 
and kinetic energy remains constant in an adiabatic flow. Enthalpy performs a 
similar role that internal energy performs in a nonflowing system. The difference 
between the two types of systems is the flow work pv required to push the fluid 
through a section. 

Bernoulli and Euler Equations For inviscid flows, the steady form of the 


momentum equation is the Euler equation, 


— + V dV = 0 


Euler Eqn 


(14.55) 


Integrating along a streamline, we get the Bernoulli’s euqation for a 
compressible flow as 


■ dp V" 

-1-= constant 

O 2 


Bernoulli's Eqn for 
Compressible (14.56) 
Flow 


For adiabatic frictionless flows the Bernoulli’s equation is identical to the 
energy equation. To appreciate this, we have to remember that this is an isentropic 
flow, so that the Tds equation is given by 
Tds = dh - vdp 


which yields 




Then the Euler equation (14.55) can also be written as 


VdV + d /7 = 0 


Needless to say that this is identical to the adiabatic form of the energy 
Eq. (14.54). The merger of the momentum and energy equation is attributed to the 
elimination of one of the flow variables due to constant entropy. 

Momentum Principle for a Control Volume For a finite control volume between 
Sections 1 and 2 (Fig. 14.3), the momentum principle is . 

p l A l -p 2 A 2 + F= mV 2 - m V, _ J 

~ (14.57) 


or 


P\A 1 —P2^2 + — Pi t 2^2 ~ P\V\ A\ 


where F is the x component of resultant force exerted on the fluid by the walls. 
The momentumprinciple, Eq. (14.57), is applicable even when there are frictional 
dissipative processes within the control volume. 


14.6 STAGNATION AND SONIC PROPERTIES 


The stagnation values are useful reference conditions in a compressible flow. 
Suppose the properties of a flow (such as T, p, p, etc.) are known at a point. The 
stagnation properties at a point are defined as those which are to be obtained if the 
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local flow were imagined to cease to zero velocity isentropically. The stagnation 
values are denoted by a subscript zero. Thus, the stangation enthalpy is defined as 

h 0 = h + — V 1 

For a perfect gas, this yields, 

c p T 0 = c p T+j-V 2 (14.58) 

which defines the stagnation temperature. It is meaningful to express the ratio of 
(T 0 /T) in the form 


T 2 c p T 2 yRT 


or 


— = 1 + 


7-1 




(14.59) 

If we know the local temperature (T) and Mach number (Ma), we can find out the 
stagnation temperature T 0 . Consequently, isentropic relations can be used to 
obtain stagnation pressure and stagnation density as 



(14.60) 

(14.61) 


In general, the stagnation properties can vary throughout the flow field. 


However, if the flow is adiabatic, then h + —is constant throughout the flow 


(Eq. 14.54). It follows that the h 0 , T 0 , and a 0 are constant throughout an adiabatic 
flow, even in the presence of friction. It is understood that all stagnation prop¬ 
erties are constant along an isentropic flow. If such a flow starts from a large 
reservoir where the fluid is practically at rest, then the properties in the reservoir 
are equal to the stagnation properties everywhere in the flow (Fig. 14.4). 



Fig. 14.4 An isentropic process starting from a reservoir 
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There is another set of conditions of comparable usefulness where the flow is 
sonic, Ma = 1.0. These sonic, or critical properties are denoted by asterisks: p , 
p ,a , and T . These properties are attained if the local fluid is imagined to expand 
or compress isentropically until it reachers Ma = 1. 

We have already discussed that the total enthalpy, hence T 0 , is conserved so 
long the process is adiabatic, irrespective of frictional effects. In contrast, the 
stagnation pressure p 0 and density p 0 decrease if there is friction. 

From Eq. (14.58), we note that 

V 2 = 2 c p (T 0 - 7) 


or 


V = 


2yR 
Y -1 



(14.62a) 


is the relationship between the fluid velocity and local temperature (7), in an 


adiabatic flow. The flow can attain a maximum velocity of 

1 

2 


v 

r max 


W 

7-1 


(14.62b) 


As it has already been stated, the unity Mach number, Ma = 1, condition is of 
special significance in compressible flow, and we can now write from 
Eq. (14.59), (14.60) and (14.61), 



(14.63a) 

(14.63b) 

(14.63c) 


For diatomic gases, like air y = 1.4, the numerical values are 


— = 0.8333, P = 0.5282, and = 0.6339 
Po P o 

The fluid velocity and acoustic speed are equal at sonic condition and is 


-p -r 

V = a 


[yRf] m 


(14.64a) 




i 

v* = 

— RT 0 
_y +! 

2 


(14.64b) 


We shall employ both stagnation conditions and critical conditions as reference 
conditions in a variety of one dimensional compressible flows. 


14.6.1 Effect of Area Variation on Flow Properties in Isentropic Flow 

In considering the effect of area variation on flow properties in isentropic flow, 
we shall concern ourselves primarily with the velocity and pressure. We shall 
determine the effect of change in area, A, on the velocity V, and the pressure p. 
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From Eq. (14.55), we can write 
= 0 


dp 

-1 - d 


( V 


v 2 y 


or dp = - pVdV 

Dividing by pV 2 , we obtain 

dp _ dV 

p v 2 V 


(14.65) 


A convenient differential form of the continuity equation can be obtained from 
Eq. (14.50) as 

d A dV dp 

V p 

Substituting from Eq. (14.65), 

d A _ dp dp 

A 


or 


dA_ 

A 


dp 


P v 2 


1 - 


dp/dp 


(14.66) 


Invoking the relation (14.47b) for isentropic process in Eq. (14.66), we get 


dA_ _ dp V 2 
^ ~ pV 2 [ a 2 


dp 

pv 2 


[1-Ma 2 ] 


(14.67) 


From Eq. (14.67), we see that for Ma< 1 an area change causes a pressure change 
of the same sign, i.e. positive cL4 means positive dp for Ma < 1. For Ma > 1, an 
area change causes a pressure change of opposite sign. 

Again, substituting from Eq. (14.65) into Eq. (14.67), we obtain 


dA 

A 


dV 


[1 — Ma 2 ] 


(14.68) 


From Eq. (14.68), we see that Ma < 1 an area change causes a velocity change of 
opposite sign, i.e. positive cL4 means negative dFfor Ma< 1. ForMa> l,an area 
change causes a velocity change of same sign. 

These results are summarized in Fig. 14.5, and the relations (14.67) and 
(14.68) lead to the following important conclusions about compressible flows: 

(i) At subsonic speeds (Ma < 1) a decrease in area increases the speed of flow. 
A subsonic nozzle should have a convergent profile and a subsonic diffuser 
should possess a divergent profile. The flow behaviour in the regime of Ma< 1 is 
therefore qualitatively the same as in incompressible flows. 

(ii) In supersonic flows (Ma >1), the effect of area changes are different. 
According to Eq. (14.68), a supersonic nozzle must be built with an increasing 
area in the flow direction. A supersonic diffuser must be a converging channel. 
Divergent nozzles are used to produce supersonic flow in missiles and launch 
vehicles. 
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Fig. 14.5 Shapes of nozzles and diffusers in subsonic and supersonic regimes 

Suppose a nozzle is used to obtain a supersonic stream staring from low speeds 
at the inlet (Fig. 14.6). Then the Mach number should increase from Ma= 0 near 
the inlet to Ma > 1 at the exit. It is clear that the nozzle must converge in the 
subsonic portion and diverge in the supersonic portion. Such a nozzle is called a 
convergent-divergent nozzle. A convergent-divergent nozzle is also called a de 
Laval nozzle, after Carl G.P. de Laval who first used such a configuration in his 
steam turbines in late nineteenth century. From Fig. 14.6 it is clear thatthe Mach 
number must be unity at the throat, where the area is neither increasing nor 
decreasing. This is consistent with Eq. (14.68) which shows that dVcan be non¬ 
zero at the throat only if Ma = 1. It also follows that the sonic velocity can be 
achieved only at the throat of a nozzle or a diffuser. 



The condition, however, does not restrict that Ma must necessarily be unity at the 
throat. According to Eq. (14.68), a situation is possible where Ma 1 at the 
throat if dV= 0 there. For an example, the flow in a convergent-divergent duct 
may be subsonic everywhere with Ma increasing in the convergent portion and 
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decreasing in the divergent portion with Ma 1 at the throat (see Fig. 14.7). The 
first part of the duct is acting as a nozzle, whereas the second part is acting as a 
diffuser. Alternatively, we may have a convergent-divergent duct in which the 
flow is supersonic everywhere with Ma decreasing in the convergent part and 
increasing in the divergent part and again Ma ^ 1 at the throat (see Fig. 14.8). 



Fig. 14.7 Convergent-divergent duct with Ma * 1 at throat 



Fig. 14.8 Convergent-divergent duct with Ma * 1 at throat 


14.6.2 Isentropic Flow in a Converging Nozzle 


Let us consider the mass flow rate of an ideal gas through a converging nozzle. If 
the flow is isentropic, we can write 


or 

or 


— = ——. a Ma [invoking Eqs (14.5) and (14.8)] 


m = pAV 
m _ p 
A RT 

— = -^—-JyRT -Ma 
A RT 


m _ p 

a Jr 


or 


2- - Ma 
R 
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or 




or 


m 

A 


T ^ 
1 o 


y~ x [TVY Po 
T 


[invoking Eq. (14.44)] 


or 


or 


m_ _ \J P 0 M a (T 0 
A 


Ir fio 




(y+i) 

2(7-1) 


?-■ Ma 


m_ _ jy Po Ma 
A 


Ir sjr 0 


1 + ^— 1 Ma 2 


(7+1) 

2 ( 7 - 1 ) 


(14.69) 


In the expression (14.69), p 0 , T 0 , y and R are constant. The discharge per unit 


in 

area — is a function of Ma only. There exists a particular value of Ma for which 

{mlA) is maximum. Differentiating with respect to Ma and equating it to zero, we 
get 


d {ml A) 
dMa 


y Po 


i 


1 + ^— 1 Ma 2 


-(7 + 1) 

l + ^Ma 2 '> 2(r - 1) 


or 


or 

or 


t R fio - 
(y+i) 

2 (y-l) 

1 


Ma 2 ( 7 + 1) = 2 + ( 7 - 1) Ma 2 


(7 + 1) 
2 ( 7 - 1 ) 


+ ,lio Ma 


R 4% 


-1 


7-1 


2Ma 


= 0 


Ma" (7 +1) 
2(1 + ^— 1 Ma 2 


= 0 


Ma = 1 

So, discharge is maximum when Ma = 1. 

We know that V = «Ma = .JyRT Ma. By logarithmic differentiation, we get 

dV 
V 

We also know that 

“i—i 


dMa J_dT 
Ma 2 T 


(14.70) 


1 + ^—^Ma 2 


By logarithmic differentiation, we get 


d7 

T 


( 7 -1) Ma" 


(14.59 repeated) 


dMa 


iZ—H Ma 2 Ma 


(14.71) 


2 
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From Eqs (14.70) and (14.71), we get 


d V 
V 


d Ma 
Ma 


{(y-l)/2}Afe 2 


1 


(7 ^ Ma 2 


dV 

V 


1 


1 + 


(7-1) 

2 


Ma 2 


d Ma 
Ma 


From Eqs (14.68) and (14.72) we get 


(14.72) 


d A 1 

A (Ma 2 - 1) 


1 


1 + 



d Ma 
Ma 


d A 

(Ma 2 - 1) 

d Ma 

A 

l + (7-D 

2 

Ma 


(14.73) 



By substituting Ma = 1 in Eq. (14.73), we get cL4 = 0 or A = constant. Some 
Ma = 1 can occur only at the throat and nowhere else, and this happens only when 
he discharge is maximum. When Ma = 1, the discharge is maximum and the 
nozzle is said to be choked. The properties at the throat are termed as critical 
properties which are already expressed through Eq. (14.63a), (14.63b) and 
(14.63c). By substituting Ma = 1 in Eq. (14.69), we get 


m 

A* 




(7 + 1 ) 


7+1 


2(7-0 


(14.74) 


(as we have earlier designated critical or sonic conditions by a superscript 
asterisk). Dividing Eq. (14.74) by Eq. (14.69) we obtain 



(14.75) 


From Eq. (14.75) we see that a choice of Ma gives a unique value of A/A*. The 
variation of A/A* with Ma is shown in Fig. 14.9. Note that the curve is double 
valued; that is, for a given value of A/A* (other than unity), there are two possible 
values of Mach number. This signifies the fact that the supersonic nozzle is 
diverging. 
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Fig. 14.9 Variation of A/A* with Ma in isentropic flow for Yg = 1.4 

14.6.3 Pressure Distribution and Choking in a 
Converging Nozzle 

Let us first consider a convergent nozzle as shown in Fig. 14.10(a). Figure 
14.10(b) shows the pressure ratio p/p 0 along the length of the nozzle. The inlet 
conditions of the gas are at the stagnation state (p 0 , T 0 ) which are constants. The 
pressure at the exit plane of the nozzle is denoted by p E and the back pressure is 
p B which can be varied by the adjustment of the valve. At the condition/^ =p E = 

p B , there shall be no flow through the nozzle. The pressure is p 0 throughout, as 
shown by condition (i) in Fig. 14.10(b). As p B is gradually reduced, the flow rate 
shall increase. The pressure will decrease in the direction of flow as shown by 
condition (ii) in Fig. 14.10(b). The exit plane pressure p E shall remain equal to 
p B so long as the maximum discharge condition is not reached. Condition (iii) in 
Fig. 14.10(b) illustrates the pressure distribution in the maximum discharge 
situation. When (mlA) attains its maximum value, given by substituting Ma = 1 
in Eq. (14.69), p E is equal to p . Since the nozzle does not have a diverging 
section, further reduction in back pressure p B will not accelerate the flow to 
supersonic condition. As a result, the exit pressure p E shall continue to remain at 
p even though p B is lowered further. The convergent-nozzle discharge against 
the variation of back pressure is shown in Fig. 14.11. As it has been pointed out 
earlier, the maximum value of (mlA) at Ma = 1 is stated as the choked flow. With 
a given nozzle, the flow rate cannot be increased further, thus neither the nozzle 
exit pressure, nor the mass flow rate are affected by lowering p B below p . 
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(b) 

Fig. 14.10 (a) Compressible flow through a converging nozzle (b) Pressure 

distribution along a converging nozzle for different values of back 
pressure 



Fig. 14.11 Mass flow rate and the variation of back pressure in a converging nozzle 

However for p B less than p*, the flow leaving the nozzle has to expand to match 
the lower back pressure as shown by condition (iv) in Fig. 14.10(b). This expan¬ 
sion process is three-dimensional and the pressure distribution cannot be pre¬ 
dicted by one-dimensional theory. Experiments reveal that a series of shocks form 
in the exit stream, resulting in an increase in entropy. 
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14.6.4 Isentropic Flow in a Converging-Diverging Nozzle 

Now consider the flow in a convergent-divergent nozzle (Fig. 14.12). The 
upstream stagnation conditions are assumed constant; the pressure in the exit 
plane of the nozzle is denoted by p E ; the nozzle discharges to the back pressure, 
p B . With the value initially closed, there is no tlow through the nozzle; the pres- 

sure is constant at p 0 . Opening the valve slightly produces the pressure 
distribution shown by curve (i). Completely subsonic flow is discerned. Then p B 
is lowered in such a way that sonic condition is reached at the throat (ii). the tlow 
rate becomes maximum for a given nozzle and the stagnation conditions. On 
further reduction of the back pressure, the flow upstream of the throat does not 
respond. However, if the back pressure is reduced further (cases (iii) and (iv)), 
the flow initially becomes supersonic in the diverging section, but then adjusts to 
the back pressure by means of a normal shock standing inside the nozzle. In such 
cases, the position of the shock moves downstream as p B is decreased, and for 
curve (iv) the normal shock stands right at the exit plane, t he flow in the entire 
divergent portion up to the exit plane is now supersonic. When the back pressure 
is reduced even further (v), there is no normal shock anywhere within the nozzle, 
and the jet pressure adjusts to p B by means of oblique shock waves outside the 
exit plane. A converging diverging nozzle is generally intended to produce 

supersonic flow near the exit plane. If the back pressure is set at (vi), the flow will 
be isentropic throughout the nozzle, and supersonic at nozzle exit. Nozzles 
operating at p B = p v , (corresponding to curve (vi) in Fig. 14.12) are said to be at 
design conditions. Rocket-propelled vehicles use converging-diverging nozzles 
to accelerate the exhaust gases to the maximum possible velocity to produce high 




Fig. 14.12 Pressure distribution along a converging-diverging nozzle for 
different values of back pressure p B 
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14.7 NORMAL SHOCKS 


Shock waves are highly localized irreversibilities in the flow. Within the distance 
of a mean free path, the flow passes from a supersonic to a subsonic state, the 
velocity decreases suddenly and the pressure rises sharply. To be more specific, a 
shock is said to have occurred if there is an abrupt reduction of velocity in the 
downstream in course of a supersonic flow in a passage or around a body. Normal 
shocks are substantially perpendicular to the flow and oblique shocks are inclined 
at other angles. Shock formation is possible for confined flows as well as for 
external flows. Normal shock and oblique shock may mutually interact to make 
another shock pattern. Different type of shocks are presented in Fig. 14.13. 





(c) Attached shock 

Fig. 14.13 



Different type of shocks 
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Fig. 14.14 One dimensional normal shock 
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Figure 14.14 shows a control surface that includes a normal shock. The fluid is 
assumed to be in thermodynamic equilibrium upstream and downstream of the 
shock, the properties of which are designated by the subscripts 1 and 2, 
respectively. 

Continuity equation can be written as 

^=p l V 1 = p 2 V 2 = G (14.76) 

A 

where G is the mass velocity kg/m 2 s. 

From momentum equation, one can write 

Pi-p 2 =r j(v 2 - V,) = p 2 V 2 2 - pi H 2 

or Pl +pi V\ =p 2 + P 2 (14.77a) 

or F x =F 2 (14.77b) 

where F = p + pV 2 can be termed as impulse function. 

The energy equation may be written as 


+ 2 ^ 2 + 2 ^ 01 ^ 02 (14.78) 

where h 0 is stagnation enthalpy. 

From the second law of thermodynamics, it may be written as 

5 2 -S!>0 (14.79) 

But Eq. (14.79) is of little help in calculating actual entropy change across the 
shock. To calculate the entropy change, we have 

Tds = dh-vdp (14.33 repeated) 

For an ideal gas we can write 


, dT dp 

ds = c„ - R—— 

p T p 


For constant specific heat, this equation can be integrated to give 


s 2 -s 1 



-R In Pl 
Pi 


(14.80) 


For an ideal gas the equation of state can be written as 

p = pRT (14.8 repeated) 

Equations (14.76), (14.77a), (14.78), (14.80) and (14.8) are the governing equa¬ 
tions for the flow of an ideal gas through normal shock. If all the properties at 
state “1” (upstream of the shock) are known, then we have six unknowns (T 2 ,p 2 , 
p 2 , V 2 , h 2 , s 2 ) in these five equations. However, we have known relationship 
between h and T [Eq. (14.17)] for an ideal gas which is given by dh = c p dT. For 
an ideal gas with constant specific heats, 
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Ah = h 2 h x =c p (T 2 - T x ) (14.81) 

Thus, we have the situation of six equations and six unknowns. 

If all the conditions at state “1” (immediately upstream of the shock) are 
known, how many possible states “2” (immediate downstream of the shock) are 
there? The mathematical answer indicates that there is a unique state “2” for a 
given state “1”. Before describing the physical picture and precise location of 
these two states let us introduce Fanno line and Rayleigh line flows. 

14.7.1 Fanno Line Flows 

If we consider a problem of frictional adiabatic flow through a duct, the governing 
Eqs (14.76), (14.78), (14.80) (14.8) and (14.81) are valid between any two points 
“1” and “2”. Equation (14.77a) requires to be modified in order to take into 
account the frictional force, R x , of the duct wall on the flow and we obtain 

R x + P\A — p 2 A = m V 2 — m V\ (14.82) 

So, for a frictional flow, we thus have the situation of six equations and seven 
unknowns. If all the conditions of “1” are known, how many possible states “2” 
are there? Mathematically, we get number of possible states “2”. With an infinite 
number of possible states “2” for a given state “1”, what do we observe if all 
possible states “2” are plotted on a T-s diagram? The locus of all possible states 
“2” reachable from state “1” is a continuous curve passing through state “1”. 
However, the question is how to determine this curve? Perhaps the simplest way 
is to assume different values of T 2 . For an assumed value ofT 2 , the corresponding 
values of all other properties at “2” and R x can be determined. 

A 


Ma < 1 



->- 

s 

Fig. 14.15 Fanno line representation of constant area adiabatic flow 

The locus of all possible downstream states is called Fanno line and is shown 
in Fig. 14.15. Point “6” corresponds to maximum entropy where the flow is sonic. 
This point splits the Fanno line into subsonic (upper) and supersonic (lower) 
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portions. If the inlet flow is supersonic and corresponds to point 1 in Fig. 14.15, 
then friction causes the downstream flow to move closer to point “b” with a 
consequent decrease of Mach number towards unity. Each point on the curve 
between point 1 and “b” corresponds to a certain duct length L. As L is made 
larger, the conditions at the exit move closer to point “ b ”. Finally, for a certain 
value of L , the flow becomes sonic. Any further increase in L is not possible 
without a drastic revision of the inlet conditions. Consider the alternative case 
where the inlet flow is subsonic, say, given the point l'in Fig. 14.15. As L 
increases, the exit conditions move closer to point “b”. If L is increased to a 
sufficiently large value, then point “b” is reached and the flow at the exit becomes 
sonic. The flow is again choked and any further increase in L is not possible 
without an adjustment of the inlet conditions. 

14.7.2 Rayleigh Line Flows 

If we consider the effects of heat transfer on a frictionless compressible flow 
through a duct, the governing Eq. (14.76), (14.77a), (14.80), (14.8) and (14.81) 
are valid between any two points “1” and “2”. Equation (14.78) requires to be 
modified in order to account for the heat transferred to the flowing fluid per unit 
mass, d Q, and we obtain 

dQ = h 02 ~^oi (14.83) 

So, for frictionless flow of an ideal gas in a constant area duct with heat transfer, 
we have again a situation of six equations and seven unknowns. If all conditions 
at state “1” are known, how many possible states “2” are there? Mathematically, 
there exists infinite number of possible states “2”. With an infinite number of 
possible states “2” for a given state “1”, what do we observe if all possible states 
“2” are plotted on a T-s diagram? The locus of all possible states “2” reachable 
from state “1” is a continuous curve passing through state “1”. Again, the 
question arises as to how to determine this curve? The simplest way to go about 
this problem is to assume different values of T 2 . For an assumed value of T 2 , the 
corresponding values of all other properties at “2” and <5 Q can be determined. 
The results of these calculations are shown on the T-s plane in Fig. 14.16. The 
curve in Fig. 14.16 is called the Rayleigh line. 

At the point of maximum temperature (point “c” in Fig. 14.16), the value of 

Mach number for an ideal gas is 1 14j ■ At the point of maximum entropy, the 
Mach number is unity. On the upper branch of the curve, the flow is always 
subsonic and it increases monotonically as we proceed to the right along the 
curve. At every point on the lower branch of the curve, the flow is supersonic, and 
it decreases monotonically as we move to the right along the curve. Irrespective 
of the initial Mach number, with heat addition, the flow state proceeds to the right 
and with heat rejection, the flow state proceeds to the left along the Rayleigh line. 
For example, let us consider a flow which is at an initial state given by 1 on the 
Rayleigh line in fig. 14.16. If heat is added to the flow, the conditions in the 
downstream region 2 will move close to point “6”. The velocity reduces due to 
increase in pressure and density, and Ma approaches unity. If 8Q is increased to a 
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sufficiently high value, then point “b” will be reached and flow in region 2 will be 
sonic. The flow is again choked, and any further increase in 8Q is not possible 
without an adjustment of the initial condition. The flow cannot become subsonic 



14.7.3 The Physical Picture of the Flow through a Normal 
Shock 

It is possible to obtain physical picture of the flow through a normal shock by 
employing some of the ideas of Fanno line and Rayleigh line Flows. Flow through 
a normal shock must satisfy Eqs (14.76), (14.77a), (14.78), (14.80), (14.8) and 
14.81). Since all the condition of state “1” are known, there is no difficulty in 
locating state “1” on T-s diagram. In order to draw a Fanno line curve through 
state “1”, we require a locus of mathematical states that satisfy Eqs (14.76), 
(14.78), (14.80), (14.8) and (14.81). The Fanno line curve does not satisfy 
Eq. (14.77a). A Rayleigh line curve through state “1” gives a locus of 
mathematical states that satisfy Eqs (14.76), (14.77a), (14.80), (14.8) and 
(14.81). The Rayleigh line does not satisfy Eq. (14.78). Both the curves on a 
same T-s diagram are shown in Fig. 14.17. As we have already pointed out, the 
normal shock should satisfy all the six equations stated above. At the same time, 
for a given state “1”, the end state “2” of the normal shock must lie on both the 
Fanno line and Rayleigh line passing through state “ 1Hence, the intersection of 
the two lines at state “2” represents the conditions downstream from the shock. In 
Fig. 14.17, the flow through the shock is indicated as transition from state “1” to 
state “2”. This is also consistent with directional principle indicated by the second 
law of thermodynamics, i.e. s 2 > .s'|. From Fig. 14.17, it is also evident that the 
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flow through a normal shock signifies a change of speed from supersonic to 
subsonic. Normal shock is possible only in a flow which is initially supersonic. 



Fig. 14.17 Intersection of Fanno line and Rayleigh line and the solution for 
normal shock condition 


14.7.4 Calculation of Flow Properties Across a Normal Shock 

The most easy way to analyze a normal shock is to consider a control surface 
around the wave as shown in Fig. 14.14. The continuity equation (14.76), the 
momentum equation (14.77) and the energy equation (14.78) have already been 
discussed earlier. The energy equation can be simplified for an ideal gas as 

7oi = T 02 (14.84) 

By making use of the equation for the speed of sound (14.49) and the equation 
of state for ideal gas (14.8), the continuity equation can be rewritten to include 
the influence of Mach number as: 

= ^-Ma 2 V7^7T (14.85) 

ivi j TVi2 

The Mach number can be introduced in momentum equation in the following way: 
Pi ^2 — Pi V\ = P\ ~Pi 


P \+ v\ =p-> + -^- L - 

RT X RT 2 
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Pi (1 + /Ma)) =p 2 (1 + yMa 2 2 ) (14.86) 

Rearranging this equation for the static pressure ratio across the shock wave, we 
get 


Pi_ = (1 + 7 Ma,) 
Pi (1 + rMa;) 


(14.87) 


As we have already seen that the Mach number of a normal shock wave is 
always greater than unity in the upstream and less than unity in the downstream, 
the static pressure always increases across the shock wave. 

The energy equation can be written in terms of the temperature and Mach 
number using the stagnation temperature relationship (14.84) as 


7 L = { l + [7-D/2]Maf 
Ti {l + [r-l)/2]Ma| 

Substituting Eqs (14.87) and (14.88) into Eq. (14.85) yields the following 
relationship for the Mach numbers upstream and downstream of a normal shock 
wave: 



Maj = 1 in Eq. (14.90) results in Ma 2 = 1. Equations (14.87) and (14.88) also 
show that there would be no pressure or temperature increase across the shock. In 
fact, the shock wave corresponding to Maj = 1 is the sound wave across which, 
by definition, pressure and temperature changes are infinitesimal. Therefore, it 
can be said that the sound wave represents a degenerated normal shock wave. 


14.7.5 Oblique Shock 

The discontinuities in supersonic flows do not always exist as normal to the flow 
direction. There are oblique shocks which are inclined with respect to the flow 
direction. Let us refer to the shock structure on an obstacle, as depicted 
qualitatively in Fig. 14.18. The segment of the shock immediately in front of the 
body behaves like a normal shock. Oblique shock is formed as a consequence of 
the bending of the shock in the free-stream direction. Sometimes in a supersonic 
flow through a duct, viscous effects cause the shock to be oblique near the walls, 
the shock being normal only in the core region. The shock is also oblique when a 
supersonic flow is made to change direction near a sharp corner. 
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Fig. 14.18 Normal and oblique shock in front of an obstacle 

The same relationships derived earlier for the normal shock are valid for the 
velocity components normal to the oblique shock. The oblique shock continues to 
bend in the downstream direction until the Mach number of the velocity 
component normal to the wave is unity. At that instant, the oblique shock 
degenerates into a so called Mach wave across which changes in flow properties 
are infinitesimal. 

Let us consider a two-dimensional oblique shock as shown in Fig. 14.19. 


Shock wave 



Fig. 14.19 Two dimensional oblique shock 

In analyzing flow through such a shock, it may be considered as a normal 
shock on which a velocity v (parallel to the shock) is superimposed. The change 
across shock front is determined in the same way as for the normal shock. The 
equations for mass, momentum and energy conservation are, respectively, 

P\ u \ = Pi u 2 (14.91) 

Pi«i («i - « 2 ) =P2~P\ (14.92) 
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Pi 


7-1 Pi 


K _ 7 Pi 

2 7 -1 Pi 


or 


7 Pi = 
7-1 Pi 2 


Pi 


(14.93) 


2 /-I p, 2 /-I p 2 

These equations are analogous to corresponding equations for normal shock. In 
addition to these, we have 

— = Ma, sin a and — = Ma-, sin B 
a x a 2 

Then modifying normal shock relations by writing Maj sin a and Ma 2 sin /) in 
place of Maj and Ma 2 , we obtain 


p 2 _ 1y Maj - sin 2 a- 7 + 1 
Pi 7 + 1 

Pi = tan /? _ 7 — 1 | _2_ 

p 2 tana 7 + 1 (7 + l)Ma 2 sin 2 a 

„2 „ ; _2 o _ 2 + (7 -1) Ma 2 sin 2 a 


Ma 2 sin" /j = 


(14.94) 

(14.95) 

(14.96) 


1 + tan^ a(tan /3/tan a) 

Note that although Ma 2 sin /3 < 1, Ma 2 may be greater than 1. So the flow behind 
an oblique shock may be supersonic although the normal component of velocity is 
subsonic. In order to obtain the angle of deflection of flow passing through an 
oblique shock, we use the relation 

n . , ^ tana- tan B 

tan 0 = tan (a- p) = 


1 + tan a tan 

tan a - (tan /3 / tan a) tan a 
1 + tan 2 a(tan /3/tana) 

Having substituted (tan /?/tan a) from Eq. (14.95), finally we get the relation 


. Ma?sin2a- 
tan 0 = 1 


2 cota 


Ma, (7 + cos2a) + 2 


(14.97) 


Sometimes, a design is done in such a way that an oblique shock is allowed 
instead of a normal shock. The losses for the case of oblique shock are much less 
than those of normal shock. This is the reason for making the nose angle of the 
fuselage of a supersonic aircraft small. 


Summary 


• Fluid density varies significantly due to a large Mach number (Ma = 
Via) flow. This leads to a situation where continuity and momentum 





























